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Summary 

Miner's  rule  for  Che  cumulative  damage  due  Co  fatigue,  a  determi¬ 
nistic  formula  which  Is  well  known  In  engineering  practice,  has  been 
examined  earlier  from  a  probabilistic  point  of  view  with  Birnbaum  in 
(2].  Here  the  assumptions  of  that  model  are  weakened.  Previously 
the  basic  assumptions  were  r.hat  crack  growth  was  stochastic  In  nature 
with  incremental  extensions  having  a  distribution  with  increasing 
failure  rate,  and  that  the  cycle  of  load  fluctuations  was  fixed  and 
then  repeated  under  program.  We  now  assume,  instead  of  the  distributions 
of  incremental  crack  extension  having  increasing  failure  rate,  only  that 
for  a  given  load  fluctuation,  the  expected  residual  damage  increment 
either  in  crack  initiation  or  extension,  given  the  damage  exceeds  a 
preassigned  amount  is  less  than  the  damage  increment  which  was  expected 
for  that  load  fluctuation  before  it  was  imposed.  We  also  weaken  the 
assumptions  concerning  the  type  of  loading  spectra  which  are  admitted, 
considering  the  case  of  random  load  fluctuations  which  are  cyclic  in 
distribution. 

Utilizing  results  from  renewal  theory  we  study  the  expected  number 
of  cycles  until  failure  under  both  programmed  and  random  loading  spectra 
and  exhibit  conditions  of  dependence  upon  load  history  under  which  a 
generalization  of  Miner's  rule  agrees  with  the  mathematical  expectation 
of  fatigue  life.  Under  other  conditions  of  dependence  we  obtain  bounds 

r  4 

for  the  expected  number  of  cycles  to  failure  under  both  programmed  and 


random  loading  spectra. 


I  V  ™ 


1 .  Introduction 

In  an  earlier  publication  [2],  Ilirnbaum  and  Saunders  gave  a  statisti¬ 
cal  interpretation  of  Miner's  rule  (a  deterministic  formula  appearing 
in  ( 6  J )  which  has  been  used  historically  under  almost  all  conditions  in 
fatigue  analysis.  In  that  paper  Miner's  rule  was  shown,  under  certain 
plausible  assumptions,  to  be  the  mathematical  expectation  of  a  stochastic 
variable  regarded  as  fatigue  life.  This  stochastic  variable  was  the 
number  of  periods  of  cyclic  oscillations  necessary  to  force  the  fatigue 
crack,  of  which  the  incremental  extensions  per  cycle  themselves  were 
regarded  as  random  variables  of  a  given  type,  to  exceed  a  critical  crack 
length.  This  critical  length  was  also  subject  to  chance  fluctuations  due 
to  various  causes. 

In  what  follows  we  will  consider  only  standardized  specimens  of  a 
material  which  are  subjected  to  fluctuating  stresses  due  to  loading  of  a 
periodic  or  stochastic  nature.  To  be  more  specific,  for  a  load  (or  load 
function)  we  mean  a  continuous  piecewise  linear  function  on  the  positive 
real  axis,  the  value  of  which  at  any  time  gives  the  stress  imposed  by  the 
deflection  of  the  material  specim3i.  This  loading  function  will  be  genericallv 
denoted  by  A,  with  or  without  affixes.  Thus  a  load  function  determines 
such  parameters  as  maximum  stress,  minimum  stress  and  average  stress 
which  are  usually  used  to  define  each  loading  oscillation. 

A  detailed  comparison  has  been  made  in  [2]  among  several  sets  of 
assumptions  and  their  consequences  for  Miner's  rule  in  its  traditional 
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form.  We  quote  here  only  rlie  most  frequently  used  deterministic  form 
of  Miner's  rule,  namely  that  fatigue  life  under  a  spectrum  of  loads  can 
be  expressed  as  a  harmonic  mean  of  the  lives  under  the  repetition  of 
certain  fixed  loads  which  comprise  the  spectrum.  This  form  is  dependent 
upon  the  assumption  that  the  order  of  the  load  oscillations  can  be  per¬ 
muted  in  any  cycle  with  the  same  resulting  fatigue  damage.  Thus  it  be¬ 
comes  necessary  to  count  only  the  number  of  oscillations  of  each  load  of 
a  given  kind. 


If  a  given  cycle  contains  various  numbers  of  different  oscillations, 
say  n^  oscillations  of  the  i1"'1  load  type  among  k  distinct  types,  then 
the  number  of  such  cycles  which  can  be  repeated  until  failure  is 


(1.1) 
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where  N  is  the  number  of  oscillations  to  failure  under  repeated 
application  of  the  same  ic"  load. 


By  this  cumulative  damage  rule  we  determine  the  fraction  of  damage 
accrued  during  one  cycle  and  use  its  reciprocal  to  calculate  the  total 
life.  In  practice  are  determined  from  available  data  on  the  regres¬ 

sion  of  stress  versus  number  of  cycles  to  failure.  The  n^  are  calculated 
from  a  typical  spectrum  of  loads  during  the  cycle.  Then  Miner's  rule  in 
the  form  given  in  (1.1)  is  used  to  determine  the  life  N  in  cycles  to 
failure . 


ft 


The  objection  that  has  been  made  to  Miner's  rule  is  that  under  certain 
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programmed  loads  witliin  laboratory  control,  the  rule  can  predict  life 
either  conservatively  or  unconservatively  depending  in  part  upon  the 
order  in  which  the  loads  are  applied.  This  is  to  say  that  the  average 
life  can  be  made  significantly  less  or  greater  than  that  predicted  by 
(1.1)  by  selecting  and  repeating  certain  sequences  of  loads  of  in¬ 
creasing  or  decreasing  magnitudes.  The  evidence  for  tills  type  of 
behavior  is  well  known,  see  [7]  and  the  references  given  there. 

What  we  do  subsequently  in  this  paper  is  to  derive  a  formula  for 
the  expected  life  which  retains  the  influence  of  load  order.  The  basic 
assumption  made  previously  concerning  the  class  of  random  variab 1 es 
governing  incremental  damage  is  relaxed  here.  At  the  same  time  the 
stochastic  nature  of  the  imposed  loads  is  taken  into  account.  More¬ 
over,  our  formula  also  reduces  to  the  classical  form  of  Miner's  rule 
(J.l)  when  the  loads  are  deterministic  and  their  order  can  be  neglected. 

2  .  A  Probabilistic  Model 

As  in  [2],  it  is  assumed  that  fatigue  failure  is  due  to  the 
initiation,  growth  and  ultimate  extension  of  a  dominant  crack.  At 
eacli  oscillation  of  the  imposed  stress,  this  crack  is  either  being 
formed  by  the  piling  up  of  atomic  dislocations  or  being  extended  by 
the  rupture  of  molecular  bonds.  In  any  case  the  damage  accumulates  bv 
an  amount  which  is  a  random  function  of  the  magnitude  of  the  imposed 
stress  and  the  geometry  of  the  specimen,  as  well  as  the  inhomogeneity 
of  the  material  and  the  influence  of  environment.  The  incremental 
damage  at  each  fluctuation  is  therefore  a  non-negative  random  variable 
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whose  distribution  may  depend  upon  several  unknown  parameters,  the 
nature  of  which  we  do  not  specify  at  present.  In  what  follows  we 
shall  speak  only  of  "crack  extension"  but  it  is  to  be  kept  in  mind 

that  we  refer  not  only  to  the  crack  growth  phase  I n  the  usual  sense  i 

I 

but  the  sub-microscopi c  phase  of  crack  initiation  as  well. 

Let  V*  denote  a  through  l  load  fluctuations. 

That  is  ■ 1  denotes  that  portion  of  the  load  function  \  which 

extends  from  time  zero  when  the  K  id  was  relaxed  and  the  stress  zero  ; 

until  the  time  when  the  i*-"  fluctuation  has  occurred. 

Our  first  assumption  is: 

1°  The  i-k  incremental  crack  extension  during  the  last  fluctuation 
of  the  loading  history  \*  is  a  non-negative  random  variable 
Zi(Xi),  depending  only  upon  A* .  The  for  i  '  1 

are  statistically  independent  random  variables. 

This  assumption  implies  the  statistical  independence  of  the  crack 
extensions  in  each  fluctuation  not  only  from  each  other  but  from  the 
total  crack  length  as  well.  Of  course  we  do  not  preclude  functional 
dependence  between  the  successive  distributions.  In  this  manner  the 
dependence  upon  the  order  of  the  loads  is  retained.  This  assumption 
appears  to  be  sufficient  in  most  cases.  Certainly  such  dependence  is 
realistic  in  the  initial  stages  of  fatigue  crack  growth,  and  mr.y  be  even 
in  situations  such  as  those  in  which  the  stress  is  near  the  ultimate  yield 
stress  of  the  metal  or  the  crack  long  relative  to  the  specimen  under  test. 
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We  now  introduce  tlio  nomenclature:  If  X  is  a  non-negative  random 
variable,  its  complementary  distribution  function  or  co-distribution  R 
is  defined  as  unity  minus  its  distribution  or 

R(t)  =  P [ X  >  t]  for  t  >  0. 

In  life  studies  the  co-distribution  is  called  either  the  reliability 
or  the  survival  distribution.  We  have  chosen  another  name  because  such 
terms  would  be  meaningless  in  our  application  here. 

We  also  make  an  assumption  about  the  probabilistic  behavior  of 
the  incremental  crack  extension  during  a  given  load  fluctuation, 
considering  that  the  crack  extension  is  influenced  by  the  loading 
history  of  preceding  fluctuations.  Specifically  we  assume 

2°  The  incremental  growth  random  variable  7.0) » 

for  any  loading  history  ,  has  a  co-distribution 
R ( •  : a  )  which  satisfies  the  inequality,  for  all 
x  0 

l  ,  (T 

(2.1)  R(x:\)  f  P.(t:\)dt  /  R(t+x:\)dt. 

'  0  0 

The  inequality  (2.1)  is  the  definition  of  a  class  of  distributions 
described  as  "new  better  than  used  in  expectation"  and  denoted  by  the 
acronymn  NBUK.  It  can  be  interpreted  as  requiring  the  expected  residual 
growth  of  the  crack  per  fluctuation,  knowing  the  crack  extension  exceeds 
x,  is  less  than  or  equal  the  expected  crack  growth  per  fluctuation  for 
any  x  0.  We  can  write  this,  for  all  x  '  0,  as, 
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(2.2)  1  f Z (  )-x  )  xj  _  F.Z() 

whirli  is  equivalent  vith  (2.1). 

This  concept  of  NBPK  was  first  introduced  by  Barlow  and  I’roschan 
in  reliability  studies  in  11]  and  was  named  and  discussed  systematically 
by  Marshall  ant!  Proschan  in  [4]. 

We  believe  that  tie  hBUE  assumption  is  realistic  since  we  can 
now  make  one  compelling  argument  which  covers  both  the  initiation  and 
crack  extension  phases  of  fatigue.  We  have  only  to  interpret  (2.2) 
in  words  and  note  its  reasonableness;  if  a  given  amount  of  damage  is 
known  to  have  occurred  as  a  result  of  a  particular  stress  fluctuation 
then  any  expected  amount  of  damage  remaining  would  be  less  than  the 
total  amount  of  damage  expected  were  that  same  fluctuation  reimposed 
under  identical  conditions. 

In  the  previous  study,  reported  in  [2],  an  argument  is  made  that 
crack  growth  by  the  successive  rupture  of  molecular  bonds  should  be  a 
random  variable  with  increasing  failure  rate  (IFR).  As  we  now  show, 
it  would  then  t  be  WBUE. 

We  can  see  easily  that  IFR  class  contains  the  NDUF  class  if  we 
consider  the  intermediate  classification  "new  better  than  used"  (NBU) 
random  variables.  In  this  case  we  must  have  for  each  load  history 

(2.3)  R(t:  )R(x:‘)  _  R(t+x:v)  for  all  t,x  0. 

By  integrating  both  sides  of  (2.3)  with  respect  to  t  we  obtain 
the  NBUE  condition  of  (2.1).  The  inequality  (2.3)  means  that  the  random 
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residual  crack  growth  knowing  the  crack  has  already  extended  a  length 
x  >  0  during  a  fluctuation  is  stochastically  smaller  than  the  random  crack 
growth  from  the  beginning  of  the  fluctuation.  The  relationship  (2.3) 
above  can  be  rewritten  as 


(2. A) 


R(x:  )-R(t+x:A) 

R(x:> )  -  1 


R(t :  A ) . 


The  left-hand  side  of  (2. A)  being  an  increasing  function  of  x  >  0 
for  each  t  >  0  is  equivalent  with  the  general  definition  of  IFR  random 
variables,  namely,  that  -;n  R(t:\)  is  convex  for  t  0. 


Remark:  It  follows  from  1°  that  the  total  crack  length  at  the  end 

of  the  n^  cycle,  say  S^,  is  a  random  variable  which  is  the  sum  of 
independent  random  variables.  Each  summand  of  represents  the  crack 

growth  during  the  corresponding  cycle. 


We  make  another  assumption  which  replaces  the  classical  one,  at 
present  considered  to  be  erroneous,  that  the  load  oscillations  may  be 
permuted  in  any  cycle  without  altering  the  resulting  fatigue  damage, 
see  [7]. 

This  assumption  allows  the  resolution  of  a  load  history  into  an 
equivalent  one  in  terms  of  fatigue  data. 

3°  There  exists  a  finite  set  of  loading  oscillations,  say 

such  that  for  any  admissible  loading 
history  there  exists  an  equivalent  ,.i  t  , ,  written 
\  for  which  in  distribution 


(2.5) 


Z(\)  =  Z(u  ). 
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In  particular  o0  is  the  loading  oscillation  such  that 

zu0)  =  0. 

When  (2.5)  holds,  the  incremental  stochastic  crack  growth 
resulting  from  the  last  oscillation  of  the  load  history  A  is  the 
same  in  distribution  as  the  incremental  crack  grot^th  following  the 
single  oscillation  u.  preceeded  by  any  number  of  repetitions  of 
itself.  Thus  ft  determines  a  partition  of  the  set  of  admissible 
loading  histories  into  equivalence  classes.  Since  ft  is  finite, 
there  exists  a  fixed  number  of  previous  oscillations  beyond  which 
the  "memory"  of  Z  does  not  extend.  Also  note  that  if  ^ 
i.e.  both  equivalent  to  the  same  element  in  ft,  it  does  not  follow 
that  the  two  histories  are  of  the  same  length. 

The  determination  of  ft  is  not  a  mathematical  problem.  Its  deter¬ 

mination  can  he  regarded  as  the  subject  of  much  of  the  recent  research  on 
the  influence  of  load  order  on  crack  growth.  In  particular  work  given  in 
[3],  [5]  and  [8]  would  be  of  that  nature. 

We  now  make  an  assumption  concerning  the  critical  crack  length,  the 
occurrence  of  which  defines  failure.  This  critical  length  can,  in  practice, 
mean  anything  from  "catastrophic  rupture  occurring"  to  "the  crack  becoming 
of  such  a  length  as  to  be  visually  inspectable" .  For  this  reason  we  take 
W,  the  critical  crack  length,  to  be  a  random  variable  such  that 

4°  The  critical  crack  size  W  is  statistically  independent 

of  the  crack  length  at  the  end  of  the  n^1  cycle  for 

eacli  n=l  ,2  , . . .  . 
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This  means  simply  that  knowing  the  length  W  at  which  the  crack 
will  become  critical,  i.e.  failure  will  take  place,  gives  no  information 
about  the  stochasLic  behavior  of  crack  growth  and  vice  versa. 

3.  Programmed  Loads 

Let  £  be  a  subset  of  the  continuous  real  valued  piecewise  linear 
functions  on  the  real  line  each  taking  the  value  zero  except  for  some 
interval  of  the  form  (0,m)  where  m  is  an  integer.  We  call  each  element 
X  c  £  a  load  spectrum  whenever  the  salient  points  of  X  can  occur  only 
at  the  integers  1,2,3,...  and  (primes  denoting  derivative) 

sgn  X'(x)  =  -sgn  X ' (x+1)  for  all  non-integral  x  >  0. 

The  least  positive  integer  m  such  that  X(x)  =  0  for  x  m  is 
called  the  length  of  the  spectrum.  Each  unit  interval  across  which  the 
load  function  is  not  zero  we  call  a  load  fluctuation.  The  length  is 
merely  the  number  of  load  fluctuations.  Two  successive  fluctuations, 
on  one  of  which  the  load  is  increasing,  is  an  oscillation.  We  make 
the  arbitrary  convention  that  the  number  of  positive  fluctuations  is 
the  number  of  oscillations. 

We  give  an  example  in  figure  1  of  two  spectra  each  of  two  oscillations. 


Note  that  the  load  oscillations  are  taken  here  to  be  of  the  same 
duration,  i.e.  the  frequency  is  the  same.  In  actual  fact  with  most 
imposed  loads  they  may  not  be,  but  because  of  the  lack  of  frequency  effect 
in  fatigue  cycling  we  shall  conveniently  make  this  assumption.  This 
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discussion  and  justification  has  been  made  before  by  Schijve  [8]. 

In  what  follows  we  will  require  the  definition  of  continuation 
of  load  function.  Consider  A1>X2  t£  with  the  length  of  A  ^ . 

We  define  the  binary  operation  A^*A2  on  ^  by* 

A1*A2(x)  =  A1(x)  +  A2(x-m1) 

if  A1(m1~l)  and  A2(l)  are  of  same  sign  and 

Ai*A2(x)  =  A-^x)  +  A 2  (x-mj+1) 

if  X1 (mj-1)  and  A2  (1)  are  of  different  sign. 

We  define  a  programmed  load  as  the  repetition  of  a  spectrum.  It 

is  thus  a  cyclic  continuous  real  valued  piecewise  linear  function  on 

the  real  line  of  infinite  length. 


Let  A  c  £  be  given  with  length  m  >  0.  Then  A_  is  the  propvammed 
load  with  spectrum  A  whenever  _A  =  A^*A2*  •••  where  A^  =  A  and 
A^+^(x)  =  A ; (x-m)  for  i  =  1,  2,  .... 


Thus 


,th 


j  represents  the  j  cycle  of  the  programmed  spectrum.  Let 

Aj+i  be  the  partial  repetition  up  through  the  itl1  load  fluctuation  of 
st  i 

the  (j+1)  cycle,  that  is  A^+1  is  the  restriction  of  A^+1  in  domain 
to  the  interval  (jm,jm+i).  Hence  A*+1  is  the  load  history  of  the 
(j+l)St  cycle  up  through  the  i1"*1  load  fluctuation  for  each  i  >  1. 


Following  the  same  general  approach  used  previously  in  [2]  we  wish 
to  establish  bounds  for  the  expected  number  of  cyclic  repetitions  of 
A  t£  which  can  be  performed  until  failure,  as  determined  by  various 
assumptions  concerning  the  stochastic  nature  of  fatigue.  We  would  then 


-11- 


like  to  compare  these  with  the  number  given  by  Miner's  Rule, 


We  make  use  of  1°  to  write  the  random  crack  extension  under  the  j 
cyclic  repetition  of  the  load  spectrum  A,  of  length  m,  as 


th 


(3.1) 


m 

Y .(')  =  ^  Z^O1)  j=l  ,2 , 
■  i  =  l  2 


t . 


where  different  affixes  i  and  j  on  the  Z  1  s  indicate  independent 
replications  of  the  corresponding  random  variables. 


Let  Ny(_)  be  the  random  number  of  such  cycles  until  failure.  It 


is  defined  by 


(3.2)  [Nyo)  =  n]  =  [Sn_2  ( ■  )  <  W,  Sn(\)  1  W] 
where 

n 

(3.3)  S  ()  =  V  Y  (•)  for  n  _  1. 

"  j-1  J 

We  now  state  a  fundamental 

Lemma  1:  if  Yj,Y9,...  are  independent  and  identically  distributed 

non-negative  random  variables  witli  mean  and  S  =  Y,  +•••+  Y  is 

n  1  n 

independent  of  W  for  all  n*l,2,...  then  the  integer  valued  random 
variable  Ny(W)  defined  by 

(3. A)  [NV(W)  =  n]  =■  [S  .  W,  S  W] 

i  n-i  n  — 

always  satisfies  the  left-hand  inequality  below 

FW  I'W 

(3.5)  —  -  1  .  EN  (W)  -  — 

i.  —  i  —  ii 

and  if  the  Y.  are  NBIJE  it  satisfies  the  right-hand  inequality  as  well. 


r 
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Proof:  If  we  let  W  =  w  with  probability  one  then  the  inequality 
(3.5)  is  well  known  to  he  true  under  the  conditions  stated,  see  [1]. 

From  the  independence  of  W  and  S^,  (3.5)  follows  by  conditional 

expectation.  A  more  complete  exposition  of  this  point  wns  made  in 

[21. || 

By  Assumption  1°  the  Z^A*)  for  a11  1  s  1****»n,»  are  independent 

and  by  2°  are  NBUE.  Thus  by  a  known  result,  see  [41  or  as  can  be  easily 
proved  from  the  definition  (2.2), the  convolution  of  such  random  variables 
is  NBUE.  Hence  the  Y^(a)  for  1  i  1  are  all  random  variables  which 
are  NBUE.  Moreover,  they  are  independent  and  identically  distributed. 
Thus  from  Lemma  1  we  have,  setting 


w  -  EW,  v(\)  *=  enY(^)  anc*  EY j  (> )  =  (\) 


the  relation 


(3.6) 


w  ,  w 

—777  -  1  <  y(A)  <  -777 
U  (A  )  —  —  M  (  W 


From  (3.6)  there  follows 


Too+T  -  vM  -ToT  • 

If  a  spectrum  \  has  m  fluctuations  then  for  i=l,...,m 

denotes  the  m  load  histories,  which  may  not  all  be  distinct  in  their 
effect  upon  crack  growth.  In  fact,  it  is  well  accepted  that  the  crack 
does  not  grow  when  the  load  fluctuation  is  decreasing,  see  [51.  Addition¬ 
ally,  we  make  the  notational  convention  that 

(3.7)  ' 1  ~  wq  for  i  >  m. 


* 
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We  may  regard  each  e  fi  for  j*l,...,k  as  a  spectrum  consist¬ 

ing  of  a  single  oscillation  per  cycle.  In  this  case  Z(tUj)  =  Y(u>j)  and 
wc  set 


(3.8)  EY(u.)=u  v  =  ENV/  .  1-1,...,  k. 

J  3  j  Y(u>  )  ’  ’ 

This  information  is  usually  obtained  from  the  regression  plot  of  cycles 
to  failure  versus  points  in  ft.  This  is  the  so-called  S-N  diagram. 


(3.9) 


Then  by  applying  the  fundamental  lemma  to  the  spectrum  m  we  have 

— - —  <  p  <  —  1=1  k 

v  +1  -  j  -  v  J 


.1 


in  A  equivalent,  with  to..  Let 


We  now  define  the  number  of  a1 
) •!  denote  the  indicator  function  taking  the  value  one  if  true  and 
zero  otherwise.  Then 

k 

(3.10)  n  (X)  =  11a1  w  }  and  zi(A1)  =  V  zi(<o.){Ai  -  u  \ 

ii1  j=0 


Bv  taking  expectations  of  (3.1)  and  using  (3.10)  we  see 


(3.11) 


O)  =  1  u  n  (A) 


j  =  l 


j  j 


and  substituting  from  (3.3)  obtain 


k  nj()l)  k  n,  (A) 

(3* 12)  w  1  -,^  +  .  <  U(A)  <_  w  v  1 -  . 

i=1  J  J-l  'j 

By  rearranging  we  obtain  bounds  on  w/u(\)  which  when  substituted 
into  (3.6)  yield  the  following: 
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Theorem  1:  For  Programmed  Spectra.  If  the  stochastic  nature  of 
incremental  crack  growth  satisfies  Assumptions  1°  and  2°,  then  each 
cyclic  repetition  of  a  programmed  spectrum  a,  results  in  random  crack 
extensions  which  are  independent  replications  of  Y(\).  It  follows  the 
number  of  times  A  can  be  repeated  until  the  crack  exceeds  a  given 
stochastic  limit,  satisfying  4°,  has  finite  expectation  EN^ (, \ •  by 
Assumption  3°  there  exists  a  sufficiently  inclusive  set  ft  of  load 
oscillations  for  which  the  information  in  (3.3)  and  (3.10)  is  provided 
from  an  S-N  diagram  for  each  u>  e  ft. 


Bounds  for  the  expectation  are  then  given  by 


(3.13) 


k 

V 


(X) 


1  -  eny(>)  - 


k 

Y 


n^(X) 


j 


j-i  V 


Note  that  these  are  the  same  hounds  given  nreviouslv  in  [2|  except 

under  weaker  conditions  on  the  stochastic  nature  of  crack  growth  and  the 
broadening  of  the  interpretation  of  n.(A). 


4.  Convex  Spaces  of  Random  Spectra 


Having  weakened  the  basic  Assumptions  1°  and  2°  by  replacing  the 
1FR  class  of  distributions  hv  the  NBUE  class  and  given  a  broader  inter¬ 
pretation  to  the  weighting  factors  n  (\ )  appearing  in  Miner's  Rule,  we 
now  wish  to  relax  the  assumption  that  the  loads  which  are  to  be  imposed 
are  cyclic  in  the  deterministic  sense  used  hithertofore .  Consequently, 
we  make  an  assumption  more,  in  keeping  with  the  stochastic  nature  of  fatigue 
cycles  as  they  appear  in  practical  situations.  By  this  we  mean  the  cycle 


only  roughly  repeats  itself  while  the  individual  fluctuations  of  load  do 
not.  We  utilize  subsequently  the  convention  that  random  variables  will 
be  denoted  by  upper  case  letters  and  the  corresponding  lower  case  letters 
the  observed  values. 

Let  us  set 

(4.1)  i->+  =  <  •  X* :  1  (x)  _  0  for  all  x  0;. 

if  \  1/+  is  not  identically  zero  then  it  must  have  even  length 
m  _  2  so  without  loss  of  generality  assume  the  first  fluctuation 
is  on  (0,1),  >  is  increasing  thereon  and  moreover  is  then  increas¬ 

ing  on  the  intervals  (2j,2j+l)  for  j=0,...,y-  1  and  decreasing 
e-  the  intervals  (2j-l,2j)  for  j*l . . 

As  a  consequence  of  (4.1),  can  be  characterized  by  the  vector 

(4.2)  (  (1),  (2) . • (m-1 ) ) 

where 

0  _  (2j)  (2.1+1)  1  (2 j+2)  j-0 . f  -  1. 

A  typical  element  in  is  graphed  below  in  figure  2. 


-16- 


If  we  denote  £  =  {A  c.  £:  -A  e  .£+} ,  then  a  space  which  represents 

the  ground-air-ground  cycle  in  certain  aeronautical  fatigue  studies  is  a 
subset  of 

(4.3)  £  *£~. 

A  typical  element  of  this  space  is  given  below  in  figure  3. 


Note  that  the  continuation  is  well  defined  for  the  stresses  of  the 
cycle  during  take-off,  flight  and  landing  for  a  particular  lower  wing 
station  of  the  airframe.  However,  insofar  as  possible  we  shall  leave 
the  exact  nature  of  £  unspecified  in  order  to  obtain  as  much  generality 
as  possible. 

Let  A  denote  a  random  function  taking  values  in  £  and  A* 
denote  the  restriction  of  A  in  domain  to  the  interval  (0,i). 


Again  we  call  A*  the  ittl  (random)  history  of  the  load  and  it 
corresponds  to  A  as  previously  defined.  The  length  of  the  cycle, 
call  it  M,  may  be  a  random  integer.  For  i  <_  M,  A1  is  well  defined 
and  for  i  M  we  follow  the  convention  of  (3.7)  and  set  A1  ~  uj  . 
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We  define  the  random  load  _A  to  be  cyclic  in  distribution  when¬ 
ever 

■  =  \  * 

—  1  2 

where  x  for  i  =  1,  2,...  are  independent  and  identically  distributed 
replications  of  some  random  load  function  A  on  £  .  Note  in  the  case 
M  =  m  with  probability  one  that  we  have  equality  in  distribution  for  all 
j  #  k  of  ...  *  for  each  i  =  l,...m.  This  replaces  the  cyclic 
assumption  of  Section  3. 


During  the  cycle  let  Zj(Aj)  be  t*ie  ran^om  microscopic  crack 
extension  due  both  to  that  portion  of  the  j*1*1  random  cycle  up  to  the 
i1"1'  fluctuation,  namely  /.*,  and  the  stochastic  variation  within  the 
material,  i.e.  its  inhomogeneity. 


We  see  the  total  damage  in  the  sense  of  crack  extension  during  the 
jCl  cycle  of  load  fluctuations,  from  our  assumption,  is 


(4.4) 


Y.(A.)  *  V  zf(A*) 
J  J  i^l  J  .  J 


We  now  let 


(4.5)  VA>  -  1  Y  (A  ) 

J-l  1  3 

ho  the  total  crack  length  at  the  end  of  n  cycles.  Note  that  Sn(A_), 
by  assumption,  is  the  sum  of  n  independent  and  identically  distributed 
random  variables. 


What  we  would  like  to  prove  is  that  (3.5)  holds  for  the  random 
variable  Ny(\)’  which  is  defined  by  the  formula  corresponding  to  (3.2) 
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wit  h  S  now  defined  by  (4.5).  However,  we 
n 

as  before  because  the  assumption  of  Z^(Xj) 
to  imply  that  Z*(.\j)  are  tllc  same  clflSS* 


conditional  expectation 


cannot  proceed  exactly 
being  NBUE  is  not  sufficient 
One  may  see  this  from  the 


P[Zj(Aj)  x]  =  ER(x:A*) 

by  knowing  mixtures  of  NBUE  variables  are  not  necessarily  NBUE.  See  the 
example  in  Appendix  A.  Moreover,  the  Zi(Ai)  need  not  even  be  independent 
since  the  A*  are  not.  Therefore,  the  hope  that  the  NBUE  property  could 
be  invoked  to  use  the  right-hand  inequality  of  (3.5)  is  vain.  Con¬ 
sequently,  we  must  impose  further  assumptions  on  the  basic  model. 


We  now  state  the  crucial  assumption 


5°  For  any  x  >  0,  a  co-distribution  R(x:*)  is  a  convex 
function  over  the  convex  sample  space  £  . 

Clearly  +  and  £  are  convex  spaces.  They  represent  loading 
functions  in  which  the  specimen  is  exclusively  in  either  compression  or 
tension.  However,  spectra  which  contain  fluctuations  which  are  both 
compressive  and  extensional  may  not,  in  general,  form  a  convex  space 
unless  other  assumptions  are  made.  Thus  to  assume  that  the  spectra  of 
the  ground-air-grounu  cycle  is  a  convex  space  it  would  be  necessary  to 
assume  that  each  of  the  three  portions  would  be  of  fixed  length.  Of 
course,  this  might  be  suitable  for  some  applications  but  not  for  others. 


Before  we  proceed  with  determining  the  implications  of  this  assump¬ 
tion  let  us  discuss  its  reasonableness.  That  £  is  a  convex  space  is 
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j 

i 


merely  a  mathematical  restriction  on  the  type  of  admissible  loading. 

It  is  made  only  to  facilitate  our  analysis.  However,  the  second  part 
of  the  assumption,  to  wit  ,  R(x:0  is  a  convex  function,  is  a  strong 

inference  about  the  nature  of  fatigue  crack  growth.  Thus  any  conclusion 
we  draw  must  reflect  any  uncertainty  that  we  have  in  this  assumption. 

Consider  the  vector  representation  of  X  as  given  in  (4.2) 

It  is  clear  that  for  any  x  >  0,  R(x:X)  should  be  convex  increasing 
as  a  function  of  each  of  the  fluctuation  peaks  (maximum  stress  per 
oscillation)  X(2i+1),  i«=0,...,y-  1  and  convex  decreasing  as  a 

function  of  the  fluctuation  troughs  (minimum  stress  per  oscillation) 

A  (2 j ) ,  j=l . ^  -  1.  See  References  [3]  and  [9]  and  Appendix  B. 

What  we  assume  in  the  second  part  of  4°  is  the  joint  convexity 
over  all  the  variables  X(i),  i»l,...,m-l  restricted  by 

X  (2j  )  <  X  (2j+l)  >  X  (2j+2) ,  j-0 . f  -  1. 

Thus  our  simply  stated  assumption  of  joint  convexity  is  not  too  great 
a  conceptual  step  from  the  known  convexity  in  each  variable  separately. 

Because  £  is  assumed  to  be  a  convex  space  it  follows,  since  A 
is  a  random  function  taking  values  in  £,  that  EA  is  a  spectrum  in 
£ .  The  other  part  of  5°,  namely  that  R(x:p)  is  convex  over  £  is 
equivalent  with  the  stochastic  inequality 

(4.6)  Y(A)  _>  Y (EA)  . 

To  see  this  note  by  the  theorem  on  conditional  probabilities  that 


for  all  y  >  0 
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(A.  7)  P  t  Y  (A )  _>  y]  -  KR(y:A)  R(y  :  EA )  -  P[Y(EA)  ^  y]. 

By  adding  independent  identically  distributed  variates  each  pair 
satisfying  (A. 6)  we  have  from  the  definition  in  (A. 7)  the  stochastic 
inequality 

S  (A)  >  S  (EA) 
n  —  —  n  — 

By  analogy  with  the  fundamental  lemma  equation  (3. A)  we  have 

<NY(A)  -  nl  - 

thus  there  follows  the  stochastic  inequality 


(A. 10) 
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Thus  we  obtain  the  expression 


U  -  EY(fl)  -  ^  '  u  p*  - 


k 

S 


ill  j-1 


JKj 


k 

-  l,l"l 

j-1  J  J 


where 


(4.11) 


n  =  p*  «  E  £  {A1  =  10  }  ■  En.  (A) 

J  JV[  J  J  J 


is  the  expected  number  of  histories  A  in  the  cycle  A  equivalent 
with  w  ,  to  use  the  notation  of  (3.10).  By  multiplying  (3.9)  by 
n^  and  summing  we  have 


k  n  k  n 

w  ^  <  u  <  w  \  . 

I  v1 "  ‘  A  vj 

Rewriting  the  above  we  obtain 


1 


*  £  < 


k  n,  —  y  ~  k  n, 

V  _L  V  — L. 

T  vj  T  V1 


By  Lemma  1  we  know  that  for  every  sequence  of  Independent  identically 
distributed  random  variables,  in  particular  Yj(Aj),  j-1.2,...,  we  have 


1~ENy(A) 


Thus  we  obtain 


V  !i 

j-1  Vj 


1  1  (A)  * 


It  is  known  that  for  some  MNBUE  variables  Y,  we  have  ENy  w/EY. 
We  must  make  use  of  the  convexity  as  given  by  (4.6). 
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Combining  the  left-hand  side  of  (3.12)  with  (4.9)  we  have 


EN, 


w 


Y(A)  -  w(EA)  -  k 

I 

1 


_1 _ 

n^EA)  * 


vj+! 


This  completes  the  result  for  random  spectra.  Note  that  if  A  is 
deterministic  this  reduces  to  the  same  bounds  which  were  given  in  (3.13). 


Theorem  2:  For  Random  Spectra  on  Convex  Spaces.  If  the  stochastic 
nature  of  incremental  crack  growth  satisfies  assumptions  1° ,  2°  and  5°,  and 
each  Independent  replication  of  a  random  spectra  A  with  convex  sample 
space  results  in  independent  crack  extensions  Y(A),  then  the  number  of 
replications  of  A  that  can  be  made  until  the  crack  exceeds  a  given 
stochastic  limit,  satisfying  4°,  has  finite  expectation  EN^^  which  is 
bounded  by 


(4.12) 


1 


k  En.(A) 

V  _J - 


-  1  <  EN. 


j-1  j 


-  Y (A )  -  k  n. (EA)  * 

y  J _ 1 

j-1  Vj+1 


where  by  Assumption  3°  for  each  uij  c  fi  we  know 

(4.13)  v  -  EN  ,  and  n  (A)-  2  {X1  *  a).}  . 

J  ^  j'  J  i>l  J 


We  can  also  obtain  the 

Corollary  2.1:  Supressing  assumption  5°,  including  the  convexity 
of  the  sample  space  of  the  random  spectra  A,  we  see  that  the  number 
of  replications  of  A  until  the  crack  exceeds  the  stochastic  limit 
satisfying  4°  has  expectation  bounded  below  by  the  left  hand  side  of 
nation  (4.12). 
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5.  Conditionally  Convex  Spaces  of  Random  Spectra 

Let  us  now  consider  the  situation  when  the  loading  spectra  may  have 
its  length,  or  a  portion,  of  random  duration.  This  case  cannot  be  treated 
by  the  assumptions  used  hithertofore .  For  example,  each  portion  of  the 
ground-air-ground  cycle  such  as  take-off,  flight  and  landing  may  well  be 
of  such  significant  variability  that  it  should  be  considered  of  random 
duration.  In  this  situation  the  sample  space  of  load  spectra  would  not 
be  convex.  Nevertheless,  there  is  a  set  of  conditions  which  if  known 
(for  example,  the  gross  take-off  weight,  barometric  pressure,  and  the 
length  of  time  of  flight)  make  all  the  random  spectra  with  these  given 
boundary  conditions  have  the  same  distribution  on  a  convex  subset  of  the 
sample  space. 

Thus  we  postulate,  there  exists  a  random  couple  (A,<J>)  having  a 
joint  distribution  on  £  '  T  for  which  the  conditional  random  function 
A  |  .}> ,  for  each  t  T,  is  a  random  spectrum  taking  values  in  a  convex 
subspace  of  J?.  The  case  we  have  in  mind  is  A 1 1>  being  a  random  spectrum 
of  fixed  length  in  each  of  its  portions  as  was  considered  in  the  preceding 
section. 


What  we  wish  to  determine  is  bounds  on 


Here  again 
immediately 


eny(a,*)  where  Y(A»$)  =  1 

Z^CA1,*)  are  the  random  incremental  crack  extensions 
following  the  i^'  fluctuation  within  the  history. 


By  assumption  Yj(Aj’^j>’ 
distributed  random  variables  wi 


j=l,2,...,  are  independent  and  identically 
th  expectation  u  ■  EY(A,iO. 


r 


Rearranging  terms  using  the  properties  of  conditional  expectation 
and  with  an  obvious  modification  of  (4.10),  we  have 


k 

=  E.  V  ui:A(..Aj ;)  -  V  u  E  \  Pi (i )  . 
‘  i_l  j-1  J  i>l  J 

We  now  have 

k  k 

i;  *  ^  u  .E.H,n  (A  ]  4>)  =  V  u  n* 
j-1  J  j-1  2  3 


where  we  set,  by  analogy  with  (3.10)  and  (4.11) 


n*  *  Enj (A  ,4)  j  -  1, . . .  ,k 


Proceeding  as  before,  from  (3.9),  we  obtain 


_1 

k  nf 

v  -L 
J-i  VJ 


-  u  ~  k  n' 

A  v,+i 

j-i  j 


By  the  fundamental  lemma,  for  any  non-negative  random  variable 


1 


k  n* 

V  -J. 
J-1  "J 


-  1  <  -  -  1  <  EN  , .  .v 

“  n  -  Y(A  ,<!) 


Consider  the  set 


9*  -  W  t  3>:  ER(-:EA|-»  >  R ( • : EA  |  4> )  > 

where  an  inequality  between  functions  indicates  the  corresponding 
inequality  between  functional  values  for  all  values  of  the  domain. 
If  3**  is  not  empty,  it  follows  by  Zorn’s  lemma,  since  linearly 
ordered  subsets  of  T* ,  ordered  by  the  function  R(*:EA|<J>)  for 

.  a 

I  t  f*,  have  lower  bounds,  that  there  are  minimal  elements  of  V 
We  then  pick  <p  as  the  maximum  likelihood  of  the  minimal  elements 
rated  by  the  marginal  density  of  :. 
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Proceeding  in  the  same  manner  as  before,  by  using  the  properties 
of  conditional  expectation  and  the  convexity  of  the  domain  as  well  as 
that  of  the  co-distribution  function,  we  have 

ER(*  :A ,<t)  >_  ER(-:EA|$)  _>  R(.;EA|$). 

Ultimately  we  obtain 


(5.1) 


k 

V 


En^A.t) 


J-l  J 


1  1  E^y(A  ,1)  1 


k 

V 

J-l 


nj (EA |  IY 

Y‘ 


Theorem  3:  For  Random  Spectra  on  Conditionally  Convex  Spaces.  If  the 
stochastic  nature  of  incremental  crack  growth  satisfies  assumptions  l6, 
2°  and  3°  and  each  independent  replication  of  the  random  couple  (A,t) 
causes  an  independent  crack  extension  Y(A,t)  and  the  conditional 
spectrum  A  ]  4>  has  a  convex  sample  space  for  each  $  and  the  marginal 
distribution  is  known,  then  the  number  of  replications  of  (A ,?)  that 
can  be  made  until  the  crack  exceeds  a  given  stochastic  limit,  satisfying 
4°,  has  finite  expectation  EN^  ^  which  is  bounded  by  (5.1)  above. 
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APPENDIX  A 

Take  R(x:')  =  exp{-x/;.0) }  for  x  >  0.  One  cheeks  easily 
that  the  NBUE  property  is  satisfied.  However,  let  A  take  two 
values,  say  and  each  with  equal  probability.  Then 

ER(x:A)  *  1i[exp(-x/'^1)  +  exp(-x/p9)  ]  . 

To  have  this  mixed  distribution  be  NBUE  is  equivalent  with 


Ii(u1+u2) 


f  exp  [~^]  +  exp  [■— M  dy 
expf-x/i^)  +  exp(-x/u9) 


for  all  x  *  0.  Upon  simplifying  we  obtain  ihe  inequality 


(P1-u2)  e_x/"2  i  (Pl-5l2)  e‘x/ul 


which  is  clearly  impossible  unless  ■  u7. 
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APPENDIX  B 


In  order  to  see  that  current  hypotheses  concerning  crack  growth 
rates  are  not  in  disagreement  with  this  assumption  we  note  that 

R(x:0  being  convex  over  £  for  each  x  >  0  implies  that 

00 

^*R(x:Odx  is  also  convex  over  £  .  This  integral  is  the  expected 
crack  extension.  In  [5]  the  "crack  growth  rate"  is  given  as  pro¬ 
portional  to 

3  1) 

f(x,y)  ■  x  (x-y)  for  x  >  y  >  0 

where  x  is  the  maximum  load  and  y  is  the  minimum  load  and  both 
are  unsealed.  We  interpret  these  two  concepts  as  being  the  same. 


The  question  now  is,  does  there  exist  a  region  of  values  of 
a,b  over  which  f  is  convex  and  if  so,  is  the  region  in  conformity 
with  those  values  obtained  by  experiment. 


A  well  known  condition  both  necessary  and  sufficient  for  the 
convexity  of  twice  differentiable  functions  is 


f  *0,  f  >  0,  ff  >  f~  , 
xx  yy  xx  yy  xy 


with  subscripts  denoting  partial  differentiation.  One  checks  that 


/  f 
xx 

/f 

/f 


2ab 

x(x-y) 


b*~-b 

(x-y)2 


-ab 

x(x-y) 
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Thus  clearly  a,b  d  1  implies  f  £  0,  f  0. 

xx  yy 

2 

We  now  check  that  f  f  -  f“  is  equivalent  with 

xx  yy  xy 

f  /{  >  a2b2/ [x*"  (b“-b)  ] .  But  note  that  xy  \ff  Ji  i  y  'J  a2 -a  +  x  *b^-b 

XX  xx 

since  1  4  b  +  a  implies  (a^-a)(b2-b)  <  a2b2.  Thus 

y  y/ a" -a  +  x^b"-b  1  yab/  \/ b“-b 

is  a  sufficient  condition  for  convexity.  However  this  last  inequality 

is  equivalent  with 

2.  >  j  b"-b 

x  ”  ab-  \/(b2-b)a2-a 

Now  1^2  0,  lienee  the  inequality  above  is  true  if  the  right  hand 
side  is  always  negative,  which  is  true  if 


But  the  square  root  of  a  number  less  than  one  exceeds  the  number. 

Hence,  a  sufficient  condition  for  convexity  is 

1  i  --  -  —  +  (l-r-)(l-^)  if  an  onlv  if  b  ^  1  +  •Ja~ . 
a  a  b  a 

Thus  if  b  *  4,  a  -  9  and  if  b  *  3,  a^4.  These  values  are 

not  in  disagreement  with  observed  values.  Also  nor*?  that  in  order  to 
obtain  simplicity  we  have  been  exceedingly  wasteful  with  our  inequality. 

For  example,  from  physical  reasons  x  must  be  bounded  and  this  restriction 
has  not  been  imposed. 
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CONCLUDING  REMARKS 

In  this  paper  we  have  made  an  attempt  to  find  physically 
realizable  conditions  under  which  both  upper  and  lower  bounds  on 
the  expected  number  of  cycles  until  failure  can  be  found  and  expressed 
by  a  formula  which  is  a  generalization  of  a  variant  of  Miner's  rule. 

We  should  not  lose  sight  of  the  fact  that  some  of  the  mathematical 
structure  imposed  was  only  necessary  to  obtain  the  upper  bound  and 
thus  insure  that  the  expected  number  of  cycles  is  given  approximately 
by  one  form  of  Miner's  rule.  Without  making  the  restrictive  assumption 
5°,  we  have  the  lower  bound  (4.11)  holding  for  arbitrary  £. 

Thus  we  conclude  that  a  variation  of  Miner's  rule  gives  a 
conservative  lower  bound  for  the  expected  number  of  cycles  until 
failure  under  very  general  conditions. 
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